# capital-coefadj.R

CoefAdj <- function(Dados, c0 = 0, cE = 1, cV = 0)
{
	E <- mean(Dados)
	V <- var(Dados)
	lab <- length(Dados)/n

	alpha <- (E^2)/V
	beta <- E/V
	c <- c0 + cE*E + cV*V

	f1 <- function (x){lab*((beta/(beta-x))^alpha)-lab-c*x}
	f2 <- function (x){(alpha*lab*((beta/(beta-x))^(alpha+1))/beta)-c}

	x1 <- x2 <- NULL
	tryCatch(
		x1 <- uniroot(f2, c(0, beta-0.00001*beta), tol = 10^-10),
		error = function(e)
			print(paste('Aviso: nao existe coef ajust para:', c0, cE, cV))
	)
	if(length(x1) == 0)
		return(NaN)
	x2 <- uniroot(f1, c(x1$root, beta-0.00001*beta), tol = 10^-10)
	x2$root
}

ProbLundberg <- function(Dados, u, c0, cE, cV) {
	r <- CoefAdj(Dados, c0, cE, cV)
	exp(-r*u)
}

